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Abstract 



> 

o 

en ■ We apply HQET to semi-leptonic B and Bg meson decays into the observed 
(N ■ 

Q\ ' charmed P wave states. In order to examine the sensitivity of the results 



to the choice of a specific model, we perform all calculations using several 
different meson models, and find that uncertainty introduced by the choice 
of a particular model is about 30%. Specifically, assuming tb = 1.50ps and 



Vcb = 0.040, we obtain branching ratios of (0.27 ± 0.08)% and (0.45 ± 0.14)% 
^ I for B DiWi and B — > -Dl^i^i decays, respectively. 

> 

X 



1 Introduction 



As more B mesons are produced at major accelerators it has become imperative to 
gain understanding of how they decay. If the final hadron state consists of a meson 
containing a heavy quark, heavy quark symmetry [jl| and heavy quark effective theory 
(HQET) 1^, 1^ provide a powerful assistance. 

Since an infinitely massive heavy quark does not recoil from the emission or ab- 
sorption of soft {E ^ Aqcd) gluons, and since magnetic interactions of such a quark 
are negligible (~ j^)^ the strong interactions of the heavy quark are independent of 
its mass and spin, and the total angular momentum j of the LDF is a good quantum 
number. Because of this, HQET leads to relations between different form factors 
describing transitions in which a hadron containing a heavy quark Q and moving 
with four-velocity v^, decays into another hadron containing a heavy quark Q', and 
moving with four- velocity v''^. In this way the number of independent form factors 
for these decays is significantly reduced. 

Semi-leptonic decays into hadrons account for over 20% of all B decays. In the 
case of B~ meson decaying into electron, neutrino, and all hadrons the branching 
ratio is 

BR{B- Xe-y) = (10.49 ± 0.46)% . (1) 

Most of the inclusive rate is accounted for hj X = D and X = D*(2010). The 
measured branching ratios for these final states are 

BR{B- ^ D*e-u) = (5.13 ± 0.84)% , (2) 
BR{B- ^ De-y) = (1.95 ± 0.55)% . (3) 

This leaves (3.4 ± 1.1)% of the hadrons unaccounted for. 

In this paper we investigate B decays into the P wave D meson states i5i(2420) 
and 1^2(2460), for which some experimental data is becoming available, and also 
the corresponding Bg decays. We use the covariant trace formalism |^, ^ and 



2 



HQET to obtain expressions for branching ratios in terms of the non-perturbative 
Isgur-Wise (IW) form factors. In order to calculate these form factors we employ 
expressions (consistent with the trace formalism), in terms of the light degrees of 
freedom (LDF) wave functions and energies . By performing all calculations using 
four different models and two different one basis state estimates, we also examine 
sensitivity of our results to the choice of a specific model. 

In Sections ^ ^ and ^ we review the covariant representation of states, decay 
rates for B —>■ Dxli^i, and the calculation of IW functions, respectively. In Section 
D we discuss the four heavy-light models which are employed in this paper. We also 
discuss estimates which use only one basis state as the wave function of the LDF. The 
model dependence of our results will be judged by the range of prediction of these 
calculations. Our main results for the S to P wave semi-leptonic branching ratios 
for B decay into Di and D2 (and corresponding Bg decays), are given in Section ^. 
Our conclusions and a comparison with experiment are summarized in Section |^. 



2 Covariant representation of states 

The covariant trace formalism, formulated in p, ^ and generalized to excited states 
in 1^, is most convenient for counting of the number of independent form factors. 
Following 1^, and using the notation of |TD|, the lowest lying mesonic states with 
mass m and four velocity v can be described as follows: 



C{v) 


1 

2^ 


/m{/) + 1)75 , 




1 

, J - 2 ' 


(4) 


C*{v,e) 


1 

- 2^ 


M(/' + i) 4 , 


J^ = l- 


1 

, J - 2 ' 


(5) 


E{v) 


1 

- 2^ 


/^i/j + l) , 




1 

, J - 2 ' 


(6) 


E*{v,e) 


1 

2^ 


/m{j6 + 1)75 4 , 


J^ = l+ 


1 

J - 2 ' 


(7) 



3 





1 

- 2^ 


3 


- 1 AY - vn] , 




= 1+ 


j 


3 

" 2 ' 


(8) 


F*{v,e) 


1 

2^ 








= 2+, 


j -- 


3 

" 2 ' 


(9) 


G{v,e) 


1 

- 2^ 








= 1" 


, 3 


3 

~ 2 ' 


(10) 


G*{v,e) 


1 

- 2^ 








= 2" 


, j 


3 

~ 2 ■ 


(11) 



In these expressions is the polarization vector for spin 1 states (satisfying e-v = 0), 



while the tensor e^'^ describes a spin 2 object (e 







0). For 



each i there are two degenerate heavy meson states (J = j ± |) forming a spin 
symmetry doublet: (C, C*) is the L = doublet, (E, E*) and (F, F*) are the two 
L = 1 doublets, and {G, G*) is an L = 2 doublet. 

In the covariant trace formalism matrix elements of bilinear currents of two heavy 
quarks (J(g) = Q'TQ) between the physical meson states are calculated by taking 
the trace (u = v ■ v'). 



{^'{v')\J{q)\^{v)) = Tt[M'{v)TM{v)]Mi{uj) 



'121 



where M' and M denote appropriate matrices from (^)-(pJ]), M = 7°M''"7°, and 
A4i{lj) represents the LDF. Again following [^, |TU[, we define the IW functions for 
the transitions of a 0~ ground state into an excited state by 



Mi(u) 



(13) 



The vector index in the last two definitions will be contracted with the one in the 
representations of excited states (PD-flTT 
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3 Decays B Dxl^i in the heavy quark limit 

Denoting the four-velocities of B and Dx mesons as and f'^, respectively, and 
assuming that lepton masses are zero, the momentum transfer is given by {uj = v-v') 

g2 = (^^^y _ rnc^v'f = {pi + p2f' = ml + m|,^ - 2mBmDx^ ■ (14) 

Using (p!^, and denoting 

X = {pi + niDxV'Y = {jubv - P2f , (15) 



the standard expression ||TT[ for the width of the semi-leptonic decay of a 5 meson 
into any of the charmed meson states Dx can be written as 

dV 



du 1287r3ml 



dx\M\^ , (16) 



where x± = mBmDxi^ i vtu^-T)- The invariant amplitude A4, 

M = ^^un''il-7')v,{Dxiv',e')\-c^,il-^')b\Biv)) , (17) 
after squaring and summing over / and Ui spins and Dx polarization, yields 

]W=IgI\V,,\'L^,H^'' . (18) 

Here, 

L'^^ = mP2-9'"'Pi-P2 + p''y2+^e''''''^PiaP2p) , (19) 
H^'' = 5:(Dx(t;',e')|c7'^(l-f )&|5(^))(^x(t^',e')|c7'^(l-7')&|5(^))^ • (20) 

pol 

The matrix elements needed in H^'^ are calculated from ( [T^ ) using (^-(|TTp, while the 
sum over Dx polarization states is performed using standard expressions for spin-1 
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and spin-2 particles, 



pol 



(21) 



pol 



\Mj^J{v)M!i\v) + lMj^J{v)Mil\v) - ImII\v)M(]^{v) . (22) 



Using kinematical identities coming from definitions (|l^) and ([T^) , and from momen- 
tum conservation, we can express in terms of u and x. Performing a simple 
integration in (p!6|), we find 



dT 



X 



Gl\V, 



where rx = mn^/mB, and the function fx is given by 

/c. (07, rcO = {to + 1)[{lo + 1)(1 - rc*f + 4^(1 - 2ujrc^ + 

(cu - l)[(cj - 1)(1 + rE*f + 4cu(l - 2ujrE* + r|0] 



[u - l){uo + l)2[(cj - 1)(1 + tfT + cu(l - 2curi. + r^)] 



(23) 



(24) 
(25) 
(26) 
(27) 
(28) 



fE'{i^,rE*) 
fFiuj,rF) 

fF*{uj,rF*) = ^{u~l){uj + iy[{uj+l){l~rF-*'f + 3uj{l~2ujrF*+r'E-.)] , (29) 

fG{uJ,rG) = '^{uj-iy{uj + l)[{uj+l){l~rGf + uj{l-2ujrG + rl)] , (30) 

/G.(^,rGO = '^{uj-iy{uj + l)[{uj-l){l + rG*)' + 3uj{l-2ujrG* + rl,)] . (31) 



Some, but not all of the above expressions can be found in earlier work ||12|-]n|. 
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4 IW functions 



The only factor in (|23|) which cannot be calculated from first principles is the IW 
function for a particular decay, . In order to estimate these form factors one has to 
rely on some model of strong interactions. In the original calculation of radiative rare 
B decays |T0[ the IW functions (0) were defined as the overlap between wave func- 



tions describing the LDF in the initial and the final mesons (AOM). These authors 
have chosen the wave functions to be eigenfunctions of orbital angular momentum L 
{a denotes all other quantum numbers), 

$aLm,(x) = R^L{r)YLmA^) , (32) 

and the form factors were given by (putting a tilde to avoid confusion with our 
definitions), 

^ci^) = {joiar))oo , (33) 

^e{uj) = y3(ji(5r))io , (34) 

U^) = v^(ji(Sr))io , (35) 

^g{oo) = V5{j2{~ar)U . (36) 



In the above 



E'V^J^ , (37) 



where = ^ denotes the "inertia parameter" of the LDF in the final heavy 
meson (with mass M' and heavy quark Q'\ The expectation values in (|33| ) - (|36|) are 
defined by 

(i^(O)La = / T''dTRl,^,{T)R^^{T)F{T) . (38) 

To calculate the above overlap integrals, in [|10[ the radial wave functions of the 
ISGW model were used. The same approach was followed in |T^ for the calcula- 



tion of semi-leptonic B meson decays into higher charmed resonances. However, by 
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comparing the covariant trace formalism of with the wave function approach 

of []TB[, it has been recently shown M that form factor definitions as pure overlap 



integrals (used in ||T0| , [T4|), are not consistent with the trace formalism. 



Under the assumption that heavy mesons can be described using simple non- 
relativistic (or semi-relativistic) quark model, the rest frame LDF wave functions 
(with angular momentum j and its projection Xj), can be written as 

^Sf W = RaL(,r)YLmA^)Xms{L,mL]^,ms\j,Xj]L,h , (39) 

where Xms represent the rest frame spinors normalized to one, xL'Xm^ = ^m'^,ms-, and 
a again represents all other quantum numbers. Performing the overlap integrals in 



the modified Breit frame (v' = — v) [|T^, using ( p9D and form factor definitions con- 
sistent with the trace formalism, one can derive expressions for the IW functions, 
their values, and values of their derivatives at the zero recoil point {u = 1), in terms 
of the LDF wave functions. Denoting [Eq = M — mq here refers to the LDF energy 
of a meson with mass M and heavy quark Q), 



(^.+^;-)>/^. (40) 



and suppressing quantum numbers a' and a, we have: 
• C ^ (C, C*) transitions. 



ec(^) = ^-(jo(ar))oo , (41) 
ic{l) = (l)oo, (42) 
ec(l) = -\-Y2^E, + E'^f<r'>oo . (43) 



C {E, E*) transitions. 



U^) = ^^ (ji(ar))io , (44) 
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c 



U^) = -(E,- + 4)(r)io, 

eUl) = -liE, + E'^{r),,-j^iE, + E:^f<r^>,, 



{F, F*) transitions. 

I 3 2 

= \l— — 7^^0iM)io 

^^(1) 



(45) 
(46) 



1 



2^3 
1 

~2VS 



2 - 1CU+ 1 

(E5 + 4)(r)io, 



(E,- + 4)(r) 
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40^3 



{Eg + E'^^<r^> 
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(47) 
(48) 
(49) 



C — > (G, G*) transitions. 

2^3 



eG(c^) 

ea(l) 



02(ar)) 



20 



20 



(50) 
(51) 



10^3 



^ {E, + E'^'{r')2o--^{Eg + EX<r'>2o • (52) 



280^3 



Note that these expressions include transitions from the ground state into radially 
excited states. If the two j = \ states are the same, E'^ = Eg and ^c(l) is normalized 
to one. Also note that from (53) and (37) it follows that the two P wave form factors 
satisfy 



and in particular 



UJ+1 

2 



V3 



(53) 



(54) 



It can be also shown that the above formulae can be generalized to any model 
involving the Dirac equation with a spherically symmetric potential. There, the wave 
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function has the form 



(55) 



where y^^j ^'^^ ^^e usual spherical spinors, A; = / (/ = j + 1) or /c = — / — 1 (/ = j — |), 



and a again denotes all other quantum numbers. Using (pSj) one finds that all the 
expressions remain unchanged, except for the expectation value (^) which 

is replaced by 

{nr))il - {F{r))^:f = I r'dr[rj;,{r)f^^{r) + g*j;,{r)g'^,{r)]F{r) . (56) 

Of course, in models with the Dirac equation the two P wave doublets are not 
degenerate any more {E'^ ^Eq), so that relations ( p^D and ( |5^ ) are no longer valid. 



5 Heavy-light models 

Although we have presented here a formalism applicable for a variety of transitions 
C C, . . . ,G*, we shall focus for the rest of the paper on the P wave transitions 
B Dilvi and B D2W1, and their B^ counterparts. The observed Di and 
D2 (or Dsi and .0*2) states are expected to be members of the j = | P wave 
doublet {F,F*), since i = \ P wave doublet {E,E*) will have 5* wave decays, and 
therefore these states should be broad resonances and correspondingly much harder 
to observe. In order to examine the sensitivity of our form factor predictions to the 
choice of a specific model and its parameters, we have performed calculations of C,f 
using four different models in the heavy-quark limit: ISGW model, semi-relativistic 
quark model (SRQM), Dirac equation with scalar confinement (DESC), and Salpeter 
equation with vector confinement (SEVC). For the ISGW model we have used the 
original parameters from |T^. For the other three models, we have fixed the string 



tension b in order to reproduce the expected linear Regge behavior, and (for a given 
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light quark mass rriu^d), varied the other parameters until a good description of 
the observed heavy-light mesons was obtained. It turns out that our results are not 
sensitive to the exact value of the light quark mass. By varying rriu^d in the range from 
0.300 GeV to 0.350 GeV we have found that specific choice of 0.300 GeV, that we 
made for the SRQM, DESC and SEVC, does not lead to significant uncertainties, as 
long as good description of the observed spin-averaged spectrum of heavy-light states 
is obtained. We have also used a single pseudo-Coulombic (PC) |T^ or harmonic 
oscillator (HO) basis wave functions to obtain form factor predictions. In this case 
the unknown energy of the LDF was estimated from the recent experimental data 
for B — » D^*Hui decays, and also from the spin-averaged masses of the known 
heavy-light states. 



5.1 ISGW model 



Because of its simplicity, the ISGW model is widely used in combination with 
HQET for calculations of different form factors. It is a non-relativistic quark model 
based on the Schrodinger equation with the usual Coulomb plus linear potential, 

V{r) = - + c + br . (57) 

3r 

The Hamiltonian of the LDF is then given by 



Hg = ^ + m^ + Vir) , (5^ 



with /i 



rriq+niQ ' 

We have used the original parameters from [p!5| , 

= O.SSGeV , 
m, = 0.55GeV , 
= l.82GeV , 
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rub = 5.12GeV , (59) 
as = 0.50 , 

c = -OMGeV , 

b = O.lSGel/^ , 

to obtain theoretical predictions for the spin-averaged heavy-hght B and D meson 
states. The resuhs are shown in Table |I]. Taking into account that the non-relativistic 
quark model with linear confinement does not yield linear Regge trajectories, this 
model provides a reasonable description of many of the known heavy-light meson 
states. Using the LDF wave functions and energies obtained from the parameters 
given above, we show in Figures |1| and |^ (with full lines) our predictions (obtained 
from d^)) for the form factor C,f in the semi-leptonic decays B Di,D2 and cor- 
responding Bg decays. For comparison we show with dashed lines the corresponding 
AOM |10|, m form factors obtained from (H). 



5.2 Semi-relativistic quark model (SRQM) 



It was observed in |TB| from Lattice QCD simulations that the ground state wave 
function describing the LDF in heavy-light mesons is in remarkably good agreement 
with the wave function that one gets from the semi-relativistic quark model. In this 
model the Hamiltonian describing the LDF is 



i/,-=^p2 + m2 + V(r) , (60) 

where 

Ary 

V(r) = — + br . (61) 

3r 

The SRQM yields linear Regge trajectories with slopes of = ^. The slope of 
the Regge trajectories in the heavy-light case is expected to be exactly twice the 
slope in the light-light case |19|, i.e. a'j^j^ = 2a^^. The observed Regge slope for 
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the light-light states is = 0.88 GeV ^ Therefore, in order to obtain the 

expected Regge behavior, we fix the string tension b to beQ 

= 0.142 GeV^ . (62) 

For a given rriu^d we vary the other parameters of the model to account for all observed 

heavy-light B and D meson states. An example of such fit is given in Table H, with 
parameters 

mu,d = OMOGeV (fixed) , 

rris = 0.512Ge\/ , 

rric = lA37GeV , 

TUb = 4.77AGeV , (63) 

as = 0.421 , 

b = O.U2GeV^ (fixed) . 

As one can see, the agreement of theoretical and experimental results is excellent. We 
show with dotted lines in Figures |^ and ^ the form factors for the decays B —>■ Di, 
and Bs Dsi,D*2, respectively. From these two figures one can see that form 
factors for the two decays are almost identical. The near equality of the B and Bg 
form factors is a refiection of the similarity of the wave functions for mesons with 
or without a strange quark. This in turns explains the near equality |Tl] of the 



rriD* — rriD — itld* — (or m^* — itlb — ttlb* — ttlb^) hyperfine differences. 



^Although this method of choosing the effective string tension ensures the correct Regge be- 
havior, it may not correspond to the correct static string tension, ff it does not, it indicates the 
interaction dynamics is incorrect. 
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5.3 Dirac equation with scalar confinement (DESC) 

Scalar confinement is the only type of confinement potential that has correct sign of 
the spin-orbit coupling. In the Dirac equation it also yields linear Regge trajectories, 
but with slope of = ^, and one can also obtain very good description of the 
spin averaged heavy-light states In this model the LDF Hamiltonian is given 



by 

H, = Ho + Pbr--^ , (64) 
Sr 

where Hq is the free particle Dirac Hamiltonian, 

Hq = a ■ p + Prriq . (65) 



Reduction of ( |64D to the set of radial equations is standard and the method of 



solution is described in |2^. In order to have the expected Regge behavior, we fix b 
to 

1 



= 0.284 GeV^ , (66) 

and, for a given rriu^d, we vary the other parameters of the model to account for all 
observed heavy-light B and D meson states. In Table ^ we show an example of such 
fit, with parameters 

mu,d = O.SOOGeV (fixed) , 
rus = 0A65GeV , 
= 1.357GeV , 
TUb = 4:.693GeV , (67) 
as = 0.462 , 
b = 0.284GeV^ (fixed) . 
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The agreement of theoretical and experimental results is again very good. The form 
factors for the decays B DijD^ and B Dsi,D*2, resulting from the DESC 
model, with parameters given above, are shown with dashed lines in Figures ^ and 
|[ Note that form factors obtained from this model are very similar to the ones 
obtained from SRQM. 



5.4 Salpeter equation with vector confinement (SEVC) 



The instantaneous version of the Bethe-Salpeter equation [|^, ^ (usually referred 
to as the Salpeter equation ||2^) is widely used for the discussion of bound state 



problems. It is also equivalent |^ to the so called "no-pair" equation |]2^, which 
was introduced in order to avoid the problem of mixing of positive and negative 
energy states that occurred in the Dirac equation for the helium atom. A similar 
problem also occurs for a single fermionic particle moving in the confining Lorentz 
vector potential. For a very long time ||28[ it has been known that there are no 
normalizable solutions to the Dirac equation in this case. 



It has been shown analytically for the heavy-light case and numerically for 
the case of fermion and antifermion with arbitrary mass E^, PU|, that in this type 
of model linear scalar confinement does not yield linear Regge trajectories. We have 
therefore used vector confinement, even though it is well known to give the wrong 
sign of the spin-orbit coupling. In terms of the free particle Dirac Hamiltonian Hq, 
potential V{r) from (^^, and the positive energy projection operator A+ defined as 



A, 



Eq + Hp 
2En 



(68) 



the LDF Hamiltonian for the heavy-light Salpeter equation with vector confinement 
is given by 



Hg = Ho + A+V{r)A+ . 



(69) 



The reduction of (|69D to a pair of coupled radial equations, as well as the solution 
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method, is described in . In Table ^ we show an example of theoretical prediction 
for the spectrum of spin-averaged heavy-light states. As in the case of SRQM and 
DESC, the agreement of theory and experiment is excellent. The parameters of the 
model were 

m„,d = 0.300Ge1/ (fixed) , 
rric = lAQAGeV , 

mfe = A.imGeV , (70) 
as = 0.534 , 
b = O.U2GeV^ (fixed) . 

Here, b was again fixed to 0.142 GeV^ since the model yields the Regge slope of 
'^'hl — ib' as in the case of SRQM. We again calculate the form factors resulting 
from the above parameters. Results for the decays B ^ Di, D2 and B —>■ Dgi, D*2, 
are shown with full lines in Figures ^ and ^ respectively. Form factors obtained from 
this model are about 10% larger than the ones obtained from the SRQM and the 
DESC. 

5.5 One basis state estimates 



Quite often (as was done in [jlO|, |3T|]) one finds in the literature estimates for form 
factors that use a single basis state as a wave function of the LDF. Usually, it is 
argued on the basis of the original ISGW model [T^, that it should be the lowest 
harmonic oscillator (HO) wave function with the scale parameter around 0.4 GeV. 



However, on the basis of lattice data [IS] one might argue that pseudo-Coulombic 



(PC) basis states |]T^ are more suitable for such a purpose. We show in Figure || a 



comparison of the lattice wave function for the heavy-light system calculated in |T^ , 
with both PC (full line) and HO IS* radial wave function with the scale parameter 
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(3s = 0.40 GeV. 

Once we choose the LDF wave function, there are still two unknown parameters 
{Eg and E'^) in the expression (^Tf) needed for the calculation of ^f- In HH the 



energy of the LDF was estimated by comparing the theoretical prediction for the 
IW function C,c obtained from (^) with the recent experimental data for the 
exclusive semi-leptonic B D'^*Hi'i decay. It is straightforward to repeat the same 
analysis using only one basis state, and obtain the range of acceptable values for 
Eg in IS* state {B,B* or D,D* mesons). The corresponding value for the other 15 
{Bg, B* mesons) or IP {Di, D2 and D^i, mesons) doublets is then determined by 
the difference between the spin averaged masses of the doublet with unknown Eq and 
the one where Eq is known {B,B* or D,D* doublets). In the following we present 
the necessary formulae for this analysis for both PC and HO wave functions. 

5.5.1 Pseudo-Coulombic basis states (PC) 

The lowest IS* and IP wave functions are ||17| 



R,s{r) = 2pfexp{-(3sr) , (71) 
Rip{r) = -^/?^/Vexp(-/3pr) . (72) 



From (|T]) we find (with Eg = E'^) 



Using this expression for with j3s = 0.40 GeV (the corresponding wave function 
is shown with full line in Figure §), and performing the analysis as described in |32], 
we find that the lowest of 0.372 per degree of freedom is obtained for 

Eg = 0.320 GeV , (74) 



17 



for the {B,B*) and {D,D*) doublets. Adding spin-averaged mass differences given 
in Tables |I|-^ we find 



Eg = 0.415 GeV , 



0.792 GeV , 
0.905 GeV , 



(75) 
(76) 
(77) 



for the {Bs,B*), {Di^D^), and {Dsi,D*2) doublets, respectively. Using these values 
for Eg, and assuming /3p = /^g = 0.40 GeV in the expressions valid for C F,F* 
transitions, 



(3T(5T{Ps + M{Eg + E'^){i^ + 1) 
((a; + ms + (3pY + (a; - l){Eg + E^f) 



eni) 



64 
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{^s + /3p)^ 

^5/2^3/2 

-IQ.^^^^XE, + E'] 



1 + 7^ 



3(i?, + 4)'' 



2 {(3s + /9p)^ 



(78) 
(79) 
(80) 



{(3s + (3pf 

we calculate form factors for B Di,D2 and B^ Dsi,D*2 decays, shown with 
the full lines in Figures ^ and [7|, respectively. 

5.5.2 Harmonic oscillator basis states (HO) 

Here, the lowest IS' and IP states. 



Risir) -- 
Rip{r) -- 

used in (^) give (with Eg = E'^] 

ec(^) = 



2f3f 
7ri/4 



exp{-f3y/2) 



n5/2 

rexp(-^|,rV2) 



3 7ri/4 



(81) 
(82) 



OJ+l 



exp 



E|(cu-1) 



^3) 
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Again, using this expression for with f3s = 0.40 GeV, and performing the analysis 
from we find that 



E, = 0.444 GeV 



(for the {B,B*) and {D,D*) doublets), yields the lowest of 0.357 per degree of 
freedom (the corresponding wave function is shown with the dashed line in Figure 
This implies 



Eg = 0.539 GeV , 



E, 



0.916 GeV , 
1.029 GeV , 



(85) 
(86) 
(87) 



for the {Bs,B*), {Di,D2), and {Dsi,D*2) doublets, respectively. The harmonic os- 
cillator expressions valid for C ^ E, F* transitions are 



^5/2^3/2 
. PP Ps (p I p/^ 



(00 + ly 



exp 



{E, + E'^\u; - 1] 
2i/3l + /3j,)iuj + 1] 



o5/2r,3/2 

I Pp Ps 

'iPl + P'pf^' 



-2 



{E, + E'] 



(E, + E': 



1 + 



1{E^ + E^ 
4 i/3l + P'p), 



(89) 
(90) 



These formulae, with f3p = Ps = 0.40 GeV and Eg values given above, yield form 
factors shown with the dashed lines in Figures ^ and |^ for B ^ Di, D2 and Bg — > 
Dsi,D*2 decays, respectively. 



6 Results for the decay rates and branching ratios 

For calculation of the decay rates we have chosen Vet = 0.040 as the reference value. 
Also, we have used fTTI r^^ = 1.50 x lO'^^^ (for B ^ Di, D* decays), and r£-^ = 
1.34 X 10~^^s (for Bs — > Dsi,D*2 decays). In order to examine sensitivity of form 
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factors to the choice of parameters of a specific model, we have fixed rric between 
1.2 GeV and 1.6 GeV, and varied other parameters of SRQM, DESC, and SEVC, 
until a good description of the spin-averaged spectra is obtained. For the one basis 



state estimates we have performed the analysis from and obtained ranges of 
acceptable Eg values {{B,B*) and {D,D*) doublets) for each Ps in the range from 
0.3 GeV to 0.5 GeV. For example, for ps = 0.3 GeV we found that Eg for {B, B*) 
and {D, D*) doublets ranges from 0.208 GeV to 0.271 GeV in the PC case, while 
the corresponding range in the HO case was 0.288 GeV to 0.374 GeV . Similarly, 
for (3s = 0.5 GeV the results were 0.346 GeV to 0.451 GeV in the PC case, and 
0.480 GeV to 0.623 GeV in the HO case. By adding the appropriate spin-averaged 
mass differences we obtained values of Eg for the other doublets. Assuming f3p = 
(3s, and varying (3s in the range from 0.3 GeV to 0.5 GeV (and Eg in the range 
corresponding to a given (3s), we obtained the acceptable ranges for all decay rates 
and branching ratios considered in this paper. 

All results obtained from different models and one basis state estimates (for C,f{^), 
^^(1), decay rates and branching ratios), are collected in Tables As one can see, 
the uncertainty introduced by the choice of parameters within a specific model is 
the largest for the SRQM (about 30%), and the smallest for the estimates which use 
pseudo-Coulombic basis states (only about 5%). For other models the uncertainty 
is about (15-25)%. We also note that the HO wave function estimates (which are 
the most commonly encountered in the literature) yield branching ratios that are 
significantly larger than those obtained from the more realistic models, for all decays 
considered in this paper. 

From the three models SRQM, DESC and SEVC, which successfully account for 
the known heavy-light masses, we obtain the following ranges for the S to P wave 
decay rates: 

2 



T{B^DiWi) = (1.16 ± 0.38) X 



V 



cb 



0.040 



lO-^'GeV^ , (91) 
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T{B ^ D;Wi) = (1.96 ± 0.66) X 
T{B, DsiWi) = (1.08 ± 0.35) x 
T{B, D*Jui) = (1.74 ± 0.62) x 



cb 



0.040 



0.040 



0.040 



Corresponding branching ratios are (using Vcb = 0.040): 



BR{B DM) 
BR{B DIW{] 
BR{Bs ^ D,M) 
BR{B, ^ D*Jui] 



(0.27 ±0.08) 
(0.45 ±0.14) 
(0.22 ±0.07) 
(0.36 ±0.13)- 



1.50ps 
1.50ps 
IMps 



% 



-% 



% 



(92) 
(93) 
(94) 

(95) 
(96) 
(97) 
(98) 



1.34ps 

For comparison with an earher work, we quote results from where B meson 
decays into charmed higher resonances were considered: 



T{B DiWi) = 0.36 X 

T{B D*Wi) = 0.52 X 
BR{B DM) = 0.0 



0.040 

V.b ' 



0.040 



IQ-^^GeV 
10-^^GeV 



BR{B D;Wi) = 0.12 



1.50ps 
1.50ps 



(99) 

(100) 
(101) 
(102) 



Let us also mention a few earlier calculations of the and form factors. In 
these form factors were computed to order 0{{uj — 1)^), 

^e{uj) = (1.43 ±0.13)- (1.86 ±0.28)(cu-l)±C((cj- 1)2) , (103) 
^piuj) = (1.14 ± 0.04) - (2.20 ±0.16)(u;- 1) ±C((cj- 1)2) , (104) 

where the quoted errors are due to a 12% variation over the scale parameter of the 
HO wave function. This should be compared with our HO estimates for ^f{^) and 
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Cf{1) given in Tables |^ and | (our errors are due to 25% variation of /^^ = /?p = 0.4). 
Assuming that the Eg value for the j = ^ P wave doublet is the same as the one 
for the j = ^ P wave doublet!, our prediction for can be obtained using (|53| ) 
and (0). In particular, using ^_f(1) = 1-23, we find ^e{^) = 1-42. Note that models 
based on the Dirac equation can result in C,e being larger than C,f- In @] the QCD 
sum rule calculation yielded ^^(l) = 1.2 ± 0.7 {C,f was not determined), while a 
Bethe-Salpeter approach of |]3^ resulted in ^£;(1) = 0.73 and ^f(I) = 0.76 

Finally, in Figures ^ and ^ we show differential branching ratios for the decays 
B DiWi and B D2W1, respectively, obtained from the four different models 
used in this paper. Differential branching ratios for the corresponding Bg decays are 
shown in Figures and |TT]. These calculations assumed the model parameters given 
in (H), (H), (H and We also assumed = 0.040, tb = 1.50 x lO^^^s and 

tb^ = 1.34 X 10~^^s. Differential branching ratios resulting from the one basis state 
estimates can be obtained using expressions (ffSf ) and (|88D, with parameters given in 

m-m and m-m. 



7 Conclusion 

In this paper we have considered semi-leptonic B and Bg meson decays into the ob- 
served charmed P wave states, in the limit where both b and c quarks are considered 
heavy. We have estimated the unknown form factors in terms of overlaps of the wave 
functions describing the final and initial states of the light degrees of freedom. Unlike 
in previous work our form factor definitions are consistent with the covariant 
trace formalism of HQET. As a result of this, we find significantly different results 

^This assumption is valid for the NRQM and the SRQM, which do not distinguish between 
doublets of the same orbital angular momentum but with different j. However, in the spirit of the 
one basis state analysis described above, where experimental mass splitting between doublets is 
used for determination of Eg, Eq's for the two P wave doublets will be slightly different. 
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for decay rates and branching ratios for processes B DiWi and B D2W1. 

In order to examine the sensitivity of our results to the choice of a specific model, 
we have performed all calculations using several different models. By fixing rric in the 
range range from 1.2 GeV to 1.6 GeV, and varying the other model parameters until 
a good description of the spin-averaged heavy-light spectrum is obtained, we have 
also examined dependence of our results on the choice of parameters within a specific 
model. We have also investigated two examples of one basis state calculations. In 
those cases the uncertainties were estimated from the range of acceptable Eg values 
consistent with the experimental data for the decay B — > D^*Hi'i (as in |^). This 
proceedure leads to the conclusion that the choice of parameters within the model 
introduces errors at the level of (5-30)%. Although six models in all are considered 
we should emphasize that three (SRQM, DESC, and SEVC), are particularly reliable 
since they account for the observed D and B spectroscopies in a very satisfactory 
manner. Between these three models, we find that the predictive accuracy for the 
unknown form factors is about 30%. 

The experimental status for B —>■ DjWi is still uncertain. At present three experi- 
mental groups have results for these decays but with possible additional non-charmed 
particle (s) X. These results are given in Table |] together with our theoretical pre- 
dictions which assume BR{Di D*7t) = 67% and BR{D* D*tx) = 20%. As 
we can see from Table ^ our predicted branching ratios are consistent with present 
measurements. In particular, all upper limits are satisfied, and where branching ra- 
tios have been determined they are somewhat greater than our predictions (in which 
there are no additional particles X). 

Finally, from Tables |^ and ^ we see that the sum of branching ratios into Di 
and D2 is about 0.72% (from the three realistic models). By counting spin states 
we estimate the total P wave meson branching ratio {E,E*,F,F*) is about 1.08%. 
Thus the P wave states account for about one third of the missing semi-leptonic B 
decays. 
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TABLES 



Table 1: Heavy-light spin averaged states. Theoretical results are obtained from the 
ISGW model |]T5[. Spin-averaged masses are calculated in the usual way, by taking 
I (|) of the triplet and ^ (|) of the singlet mass for the S{P) waves). 



State 






Spin-averaged 


Q. 


n. 


Theory 


Error 




-r D 


2^+1 L J mass (MeV) 


j 


A; 


(MeV) 


(MeV) 


cu, cd quarks 
















D(1867) C 
(2009) C 


0" 
1~ 




1^ (1974) 


1 

2 


-1 


1931 


-43 


A (2425) F 
D*{2459) F* 


1+ 

2+ 




IP (2446) 


3 
2 


-2 


2447 


1 


cs quarks 
















D,(1969) C 
D,*(2112) C* 


0" 
1- 


%" 


IS (2076) 


1 
2 


-1 


1982 


-94 


D,i(2535) F 
D*^{2573) F* 


1+ 
2+ 


'P2 


IP (2559) 


3 
2 


-2 


2473 


-86 


bu, bd quarks 
















fi(5279) C 
5*(5325) C* 


0- 

1- 




IS (5314) 


1 
2 


-1 


5188 


-126 


hs quarks 
















fi,(5374) C 
5^(5421) C* 


0" 

1- 




IS (5409) 


1 
2 


-1 


5216 


-193 
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Table 2: Heavy-light spin averaged states. Theoretical results are obtained from 
the SRQM. Parameters of the model are given in (^). Spin-averaged masses are 
calculated in the usual way, by taking | (|) of the triplet and | (|) of the singlet 
mass for the S{P) waves). 



State 


Spectroscopic label Spin-averaged 


Q. 


n. 


Theory 


Error 






mass (MeV) 


3 


k 


(MeV) 


(MeV) 


cu, cd quarks 
















D(1867) C 
D*(2009) C* 


0- 

1- 




IS (1974) 


1 

2 


-1 


1974 





A (2425) F 
D*{2459) F* 


1+ 

2+ 




IP (2446) 


3 
2 


-2 


2449 


3 


cs quarks 
















D,(1969) C 
D,*(2112) C* 


0" 
1- 




15 (2076) 


1 

2 


-1 


2075 


-1 


D,i(2535) F 
D*^{2573) F* 


1+ 
2+ 


'P2 


IP (2559) 


3 
2 


-2 


2557 


-2 


bu, hd quarks 
















fi(5279) C 
5*(5325) C* 


0- 

1" 


'So' 


IS (5314) 


1 
2 


-1 


5311 


-3 


hs quarks 
















fi,(5374) C 
S;(5421) C* 


0" 

1- 


^^1. 


15 (5409) 


1 
2 


-1 


5412 


3 
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Table 3: Heavy-light spin averaged states. Theoretical results are obtained from 
the Dirac equation with scalar confinement. Parameters of the model are given in 
(|67D . Spin-averaged masses are calculated in the usual way, by taking | (|) of the 
triplet and | (|) of the singlet mass for the S{P) waves). 



State 






Spin-averaged 


Q. 


n. 


Theory 


Error 






2^+1 L J mass (MeV) 


j 


k 


(MeV) 


(MeV) 


cu, cd quarks 
















D(1867) C 
D*(2009) C* 


0" 
1- 




IS (1974) 


1 

2 


-1 


1977 


3 


A (2425) F 
D*{2459) F* 


1+ 
2+ 


'P2 


IP (2446) 


3 
2 


-2 


2444 


-2 


cs quarks 
















D,(1969) C 
L',*(2112) C* 


0- 

1- 




IS (2076) 


1 
2 


-1 


2074 


-2 


D,i(2535) F 
D*^{2573) F* 


1+ 
2+ 


'Pl/'PI ' 

'P2_ 


IP (2559) 


3 
2 


-2 


2560 


1 


bu, hd quarks 
















5(5279) C 
5*(5325) C* 


0" 
1- 


%" 


1^ (5314) 


1 
2 


-1 


5313 


-1 


hs quarks 
















5,(5374) C 
5* (5421) C* 


0- 
1- 


%" 


1^ (5409) 


1 
2 


-1 


5410 


1 
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Table 4: Heavy-light spin averaged states. Theoretical results are obtained from the 
Salpeter equation with vector confinement (in the heavy-light limit). Parameters of 
the model are given in (^). Spin-averaged masses are calculated in the usual way, 
by taking | (|) of the triplet and j (|) of the singlet mass for the S{P) waves). 



State 






Spin-averaged 


Q. 


n. 


Theory 


Error 






2^+1 L J mass (MeV) 


j 


k 


(MeV) 


(MeV) 


cu, cd quarks 
















D(1867) C 
D*(2009) C* 


0" 
1- 




IS (1974) 


1 

2 


-1 


1980 


6 


A (2425) F 
D*{2459) F* 


1+ 
2+ 


'P2 


IP (2446) 


3 
2 


-2 


2439 


-7 


cs quarks 
















D,(1969) C 
L',*(2112) C* 


0- 

1- 




IS (2076) 


1 
2 


-1 


2072 


-4 


D,i(2535) F 
D*^{2573) F* 


1+ 
2+ 


'Pl/'PI ' 

'P2_ 


IP (2559) 


3 
2 


-2 


2564 


5 


bu, hd quarks 
















5(5279) C 
5*(5325) C* 


0" 
1- 


%" 


1^ (5314) 


1 
2 


-1 


5316 


2 


hs quarks 
















5,(5374) C 
5* (5421) C* 


0- 
1- 


%" 


1^ (5409) 


1 
2 


-1 


5407 


-2 
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Table 5: Results for the decay B Dilui obtained from four different models and 
two one basis state estimates. Errors for the SRQM, DESC, and SEVC, are due to 
variation of in the range from 1.2 GeV to 1.6 GeV . Errors for the PC and HO 
estimates are due to variation of (3s = Pp in the range from 0.3 GeV to 0.5 GeV. 
For the reference values of the B meson lifetime we take [O] t^^ = 1.50 x lO'^^^g^ 
and for the reference value of Vcb we take 0.040. 



Model 






Vet 
0.040 


' 10-^^GeV] 


BR 1-^%] 


ISGW 


0.60 


-0.89 




0.73 


0.17 


SRQM 


0.84 ±0.19 


-2.04 ±0.77 


1.01 ±0.34 


0.23 ±0.08 


DESC 


0.79 ±0.14 


-1.45 ±0.44 


1.07 ±0.24 


0.25 ±0.06 


SEVC 


1.18 ±0.19 


-4.04 ± 1.31 


1.41 ±0.22 


0.32 ±0.05 


PC 


1.42 ±0.26 


-6.23 ±2.70 


1.61 ±0.06 


0.37 ±0.02 


HO 


1.23 ±0.21 


-3.23 ± 1.18 


1.94 ±0.28 


0.44 ±0.06 



Table 6: Results for the decay B 02^1 obtained from four different models and 
two one basis state estimates. Errors for the SRQM, DESC, and SEVC, are due to 
variation of in the range from 1.2 GeV to 1.6 GeV. Errors for the PC and HO 
estimates are due to variation of Ps = Pp in the range from 0.3 GeV to 0.5 GeV. 
We take lTl| r^j^-^ = 1.50 x lO^^^^ ^^^^ y^^ ^ o.040. 



Model 




eMi) r [ 


0.040 


' 10-^^GeV] 


BR [^%] 


ISGW 


0.60 


-0.89 




1.14 


0.26 


SRQM 


0.84 ±0.19 


-2.04 ±0.77 


1.69 ±0.56 


0.38 ±0.13 


DESC 


0.79 ±0.14 


-1.45 ±0.44 


1.77 ±0.45 


0.41 ±0.11 


SEVC 


1.18 ±0.19 


-4.04 ± 1.31 


2.41 ±0.41 


0.55 ±0.10 


PC 


1.42 ±0.26 


-6.23 ±2.70 


2.80 ±0.18 


0.64 ±0.04 


HO 


1.23 ±0.21 


-3.23 ± 1.18 


3.23 ±0.55 


0.74 ±0.13 
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Table 7: Results for the decay Bs DsiWi obtained from four different models and 
two one basis state estimates. Errors for the SRQM, DESC, and SEVC, are due to 
variation of rric in the range from 1.2 GeV to 1.6 GeV. Errors for the PC and HO 
estimates are due to variation of jSs = I3p in the range from 0.3 GeV to 0.5 GeV . 
1.34 X 10-^2^ and Vcb = 0.040. 



We take 11 



ref 



Model 






0.040 


' 10-^^GeV] 


BR [%%] 


ISGW 


0.51 


-0.71 




0.53 


0.11 


SRQM 


0.84 ±0.21 


-2.09 ±0.81 


0.9f 


S±0.34 


0.20 ±0.07 


DESC 


0.82 ±0.14 


-1.59 ±0.47 


1.10 ±0.27 


0.23 ±0.06 


SEVC 


0.97 ±0.22 


-2.75 ± 1.04 


1.15 ±0.34 


0.24 ±0.07 


PC 


1.70 ±0.34 


-9.98 ±4.73 


1.57±0.11 


0.32 ±0.02 


HO 


1.42 ±0.26 


-4.59 ± 1.87 


2.01 


B±0.21 


0.43 ±0.05 



Table 8: Results for the decay B D*2liyi obtained from four different models and 
two one basis state estimates. Errors for the SRQM, DESC, and SEVC, are due to 
variation of in the range from 1.2 GeV to 1.6 GeV. Errors for the PC and HO 
estimates are due to variation of Ps = Pp in the range from 0.3 GeV to 0.5 GeV. 
For the reference values of the We take fill r]j^/ = 1.34 x 10"^^s and Vcb = 0.040. 



Model 




W) r [ 


Vet 
0.040 


' 10-^^GeV] 


BR [%%] 


ISGW 


0.51 


-0.71 




0.81 


0.17 


SRQM 


0.84 ±0.21 


-2.09 ±0.81 


1.5J 


B±0.57 


0.32 ±0.12 


DESC 


0.82 ±0.14 


-1.59 ±0.47 


1.75 ±0.45 


0.36 ±0.09 


SEVC 


0.97 ±0.22 


-2.75 ± 1.04 


1.90 ±0.60 


0.39 ±0.12 


PC 


1.70 ±0.34 


-9.98 ±4.73 


2.90 ±0.08 


0.59 ±0.02 


HO 


1.42 ±0.26 


-4.59 ± 1.87 


3.52 ±0.48 


0.72 ±0.10 
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Table 9: A summary of experimental results for B —>■ DjXe~u, where X is a possible 
non-charmed hadron. Our theoretical predictions are obtained from (^) and (|96| ) 
for the pure semi-leptonic decays B Dje~u, assuming BR{D\ D*^7i^) = 67% 
and BR{Df D*+7i-) = 20%. 



Decay Mode 


CLEO |! 


j] ALEPH 


36|] OPAL 13 This work 




[%] 


[%] 


[%] [%] 



BR{B ^ DlXe-v) < 0.67 0.51 ± 0.17 1.36 ±0.46 0.18 ±0.05 

xBR{D\ D*+7i-) (90% c.l.) 

BR{B ^ Dl^Xe'u) < 0.79 < 0.20 0.18 ±0.08 0.09 ± 0.03 

xBR{Df ^ D*+7i-) (90% c.l.) (95% c.l.) 
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FIGURES 



Figure 1: for the semi-leptonic decays B Di,D2, obtained from the ISGW 
model [|15[. The full line shows our prediction (VO), obtained from (0), while the 
dashed line is AOM prediction obtained from (135|), which is used in p^ . 



Figure 2: for the semi-leptonic decays Bg — > Dsi,D*2, obtained from the ISGW 
model |r^|. The full line shows our prediction (VO), obtained from (0), while the 
dashed line is AOM prediction obtained from (|35|) , which is used in [0, 



Figure 3: for the semi-leptonic decays B — > Di,D2, obtained from the SRQM 
(dotted line), DESC (dashed line), and SEVC (full line). Model parameters are given 
in (|63D, (|67D and ([70|), respectively. 



Figure 4: for the semi-leptonic decays Bs Dsi,D*2, obtained from the SRQM 
(dotted line), DESC (dashed line), and SEVC (full line). Model parameters are given 
in (|63D, (|67D and ([70|) , respectively. 



Figure 5: Comparison of the lattice data with the IS" pseudo-Coulombic (full hne) 
and harmonic oscillator (dashed line) wave function. For both wave functions we 
used /3s = 0.40 GeV. 



Figure 6: ^p for the semi-leptonic decays B Di,D2, obtained from one PC basis 
state (full line), and from one HO state (dashed line). In both cases we used Ps = 
Pp = 0.40 GeV. In the PC case we used Eg = 0.320 GeV {{B,B*) doublet) and 
Eq = 0.792 GeV ((£>!, /^a) doublet), and in the HO case Eg = 0.444 GeV {{B,B*) 
doublet) and Eg = 0.916 GeV {(Di^D*) doublet). 
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Figure 7: for the semi-leptonic decays Bs —>■ Dsi,D*2, obtained from one PC 
basis state (full line), and from one HO state (dashed line). In both cases we used 
= Pp = 0.40 GeV. In the PC case we used Eg = 0.415 GeV ((5„5;) doublet) 
and Eg = 0.905 GeV {{DsuD*^) doublet), and in the HO case Eg = 0.539 GeV 
{{Bs,B*) doublet) and Eg = 1.029 GeV {{DsuD*^) doublet). 

Figure 8: Differential branching ratio for the process B Dilui {C — * F) 
obtained from the four different models, with parameters given in (^) (ISGW), (|63| ) 
(SRQM), (13) (DESC), and (|70|) (SEVC). For this calculation we used V^b = 0.040 
and tb = 1.50 x 10~^^s. The kinematic limit for this decay is ujmax = 1.318. 

Figure 9: Differential branching ratio for the process B D^^li^i (C F*) 
obtained from the four different models, with parameters given in (^) (ISGW), (|63| ) 
(SRQM), (13) (DESC), and (|7D|) (SEVC). For this calculation we used = 0.040 
and Tb = 1.50 x lO^^^s. The kinematic limit for this decay is uj^ax = 1.306. 

Figure 10: Differential branching ratio for the process B^ DgiWi (C F) 
obtained from the four different models, with parameters given in (^) (ISGW), ( p3| ) 
(SRQM), (13) (DESC), and (|7|) (SEVC). For this calculation we used V^b = 0.040 
and tb^ = 1.34 x lO^^^s. The kinematic limit for this decay is uJmax = 1.296. 

Figure 11: Differential branching ratio for the process Bg -0*2^^'; (C F*) 
obtained from the four different models, with parameters given in (^) (ISGW), ( [53| ) 
(SRQM), (|63) (DESC), and (|70|) (SEVC). For this calculation we used Kfe = 0.040 
and tb^ = 1.34 x lO^^^s. The kinematic limit for this decay is Umax = 1-284. 
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Figure 9 

45 



0.025 



0.02 - 



0.015 - 



3 



0.01 - 



0.005 - 








1.05 



1.1 



1.15 
CO 



1.2 



1.25 



Figure 10 
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